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LARGE  DEFLECTIONS  OF  PLATES  WITH 
NONLINEAR  ELASTICITY  AND  HOLLOW 
SHELLS  WITH  HINGE-SUPPORTED  EDGES 


M.  S.  Kornishin,  N.  N.  Stolyarov, 
and  N.  I.  Dedov 

In  this  work,  method  of  finite  differences  is  used  to  solve 
problems  dealing  with  large  deflections  of  rectangular  plates  in 
plan  with  nonlinear  elasticity  and  hollow  shells  with  hinge- 
supported  edges,  which  are  under  the  effect  of  external  normal 
pressure  uniformly  distributed  along  the  entire  surface  or  along 
the  central  rectangular  area.  The  curvature  parameter  values 
at  which  a  crack  in  a  cylindrical  panel  occurs  are  established. 
Critical  loads  and  their  corresponding  deflections  have  been 
found  for  panels  with  different  curvature  parameters. 

1.  Principal  dependences.  Let  us  consider  a  hollow  shell 
which  is  rectangular  in  plan  with  sides  2a,  2b  and  thickness  h. 
We  will  place  the  origin  of  the  coordinates  in  the  center  of  the 
shell  and  direct  axes  x,  y  parallel  to  the  sides  of  the  plan. 

We  introduce  the  following  designations:  -  stress 

function;  u,  v,  w,  -  corresponding  displacements  of  a  point  of 
the  middle  surface  along  axes  x,  y,  z;  k^,  k 2  -  shell  curvatures 
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written  here  in  total  derivatives: 


■Sr[/f/*  -4r~|  '  di> 

-  Tz  J***»*/»{->  -/**/«  ?•(*> - 0.  ( 12 ) 

^rr  [/;/*-^r]  - p'n<*Mx)  ^  ( i3) 

The  conditions  of  articulation  (z  ■  0;  x  =  0): 


1°  /,<0)  -/♦(«). 

2°  — 'j~  0  —  condition  of  syrr.r? 


3°  —  <?*(0)  — 


•t.um 

~dx  • 


“try  of  vibration  shapes, 


Boundary  equations  when  j  /*: 


T  111 » 


lls* 


II 


o  -  absence  of  shearing  force, 

absence  of  bending  moment, 
absence  of  torque. 


Boundary  conditions  when  x  -/♦: 
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Assuming  that  =0,  from  relationships  (1)  we  obtain 
Ef(vs)  . 

f..— vt-.j  r.  — 

2(>^. 


E  )  -,-r.  F  _ 

1  11  i'/’  ~  ‘  f‘  2(i+ii) 


J  «"'£•  ~  '' 

/tVic  fi-/V 

Here  and  subsequently  1,  2  -  index  transposition  symbol 


(2) 


.0-  J 


The  expressions  for  stresses  and  moments 


m  i  a 


after  certain  transformations  with  the  use  of  relationships  (2) 
will  assume  the  form 


TV.  -  7~:  h  ,  A'  •  ~K  {■> 


(3) 


MU  -  JAJZjT)  (*/t  +JJKu)+&l%  ■  (W-  'drill'll A  +&K:, 


where 


Eh 


—  ■  * 
F/i‘ 


Eh 3 


OJ) 
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Coefficients  P,  Q,  R,  X,  Y,  Z  represent  integrals 

[Translator's  note:  equation  (5)  is  not  indicated  -  assumed  to 
be  one  of  the  following] 


P.Sh 


P - 


-O.Sn 


Jr  , 


Hi 


R= 


■ZS\ 


where 


P  - 

'  (i-f)d-v) 


*  (t-JJj(f-s) 


Having  expressed  T ^ n  terms  the  stress 

function  f , 

Pfl  r~  ^yy  <  'it  ~  @cx  1  7/z  ~  ~  ^ry  . 

from  relationships  (3)  we  obtain 


(6) 


Using  the  well  known  nonlinear  equations  of  equilibrium 
and  the  condition  of  compatibility  of  deformations  of  the 
theory  of  hollow  shells,  and  introducing  dimensionless  variables 
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,  il 


x  a  '  '  A  a  '  u  ~ k  ’  A'  Pi*  ’  f'2  hk  ' 

sts,  $  n  _  22  T  -  Tut2  T  „  ~:l  t>2 

®=fh3'  P  **~r  '  - 


r  u  __  //»£*  a/  ..  /*>/  /./  _ 

^-7^'  ' 


we  will  obtain  a  system  of  two  nonlinear  differencial  equations 

*%IIX  *  Pyyyy  '>  0.25 25 ~ 

~  ^  ^xy  ~  ^  7//,yy  ^  ^7,.ru  ^7,  *7  ^ 

^  4  4x,  ~ 2 (*+/■) \r2  ,  ^  7  ^ 

(5,.  =  p/tt  %  Vv*'„  <JX  -^x  ^ 

+A  AMft'Xx  +2XaN,2  xv  + &  f'h'Ly  • 

Presented  below  are  solution  results  of  the  geometrically  and 
physically  nonlinear  problems  of  bending  a  square  plate  and  a 
square  cylindrical  panel  using  external  normal  pressure  with 
boundary  conditions  of  hinged  support,  which  can  be  presented  as 

with  x  -  ± j 

O'  r*s  _  . 

V  "  v;T,  -a, 

with  y=±l 

<P  ~(p  =c 

rxx  yx</  l-x 


* 


Conditions  (8'-  indicate  the  equality  to  zero  on  the  con¬ 
tour  of  normal  and  tangential  stresses,  and  also  the  deflection 
and  moment . 


2.  Method  for  solving  nonlinear  problems  using  the  elec¬ 
tronic  digital  computer  and  the  calculation  results.  To  resolve 
the  problem,  we  will  use  a  method  of  finite  differences.  Having 
selected  a  rectangular  net  10  *  10,  we  substitute  the  initial 
system  of  differential  equations  (7)  with  a  system  of  differ¬ 
ence  equations,  approximating  the  biharmonic  operator  and  all 
the  derivatives  with  the  symmetric  difference  equations  with 
an  error  on  the  cr-der  of  0(h  )  where  h  -  net  spacing.  The  de¬ 
rivatives  of  function  in  boundary  conditions  we  approximate 
with  an  error  on  the  order  of  0(h  ) . 

In  accordance  with  the  boundary  conditions  the  contour 
values  of  functions  * \-0  ,  (pK=0  ,  and  the  values  of  functions 
u\.t  ,  ,  beyond  the  contour,  are  expressed  in  terms  of  the 

intracontour  values  according  to  formulas  [2]: 


with  i"  i  i 

, „  _  1 ...  _  i  a,  +±IC  ,  +  , 

*  ~  Tf  u*->  11  *•*  II  *’J  11  J 

~  3  P\-l  v  > 

with  y=±f 

= -/f  u*->  -ft  ■ 


Reproduced  from 
best  available  copy. 


Due  to  the  symmetry  of  deformational  state  of  the  shell 
relative  to  axes  x  and  y  the  number  of  difference  equations  is 
reduced  to  50.  The  obtained  system  of  50  nonlinear  algebraic 
equations  was  solved  using  a  method  of  total  iteration,  simi¬ 
larly  to  those  described  in  monograph  [2]. 


Part  of  the  calculatjon  results  with  y  =  0.3,  A  =  0, 

18-10^  for  -d'‘.  =?.  /•.  »--v7~;'  is  presented  in  the  table  and 
on  figures.  The  calculations  were  carried  out  for  the  following 
loads:  load  of  constant  intensity  distributed  throughout  the 

surface  s^  =  4ab;  with  intensity  -  along  the  cylindrical 
area  s2  =  4ab;  with  intensity  P^  along  area  =  0.36ab;  with 
intensity  along  area  s^  =  Q.04ab.  We  used  che  following 
designations  in  the  figures: 


where  f,  -  load  intensity  parameter,  -  parameter  of  total 
Joaa,  -  deflection  parameter  at  the  center. 
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The  ta’oje  shows  the  obtained,  with  the  consideration  of 
(/  -  t-y’A  ,  r  )  and  without  the  consideration  of  (}'-■)  of  physical 
nonlinearity,  parameter  values  for  upper  and  lower  (l“) 

critical  loads  and  their  corresponding  deflection  parameters  in 


FTD-HT-2 3-18 14-72 


7 


I 


center  (urf ,  t?*)  for  a  square  cylindrical  panel  uniformly  loaded 
along  the  entire  surface  in  the  plan  with  different  curvature 
parameters  and  with  ratio  A/i=0,005  •  Given  here  also  are  the 
values  expressed  in  per  cent  for  relations 


A  _J 

_  Pi, 1x1  Pi.!*! 


wo'/,  .  t= 


\ 


? 

f 


It  is  evident  from  the  table  that  physical  nonlinearity 
lowers  the  critical  loads  significantly.  With  an  increase  in 
parameter  K 2  the  effect  of  physical  nonlinearity  noticeably 
increases,  moreover,  its  effect  on  p*  is  somewhat  greater  than 
on  pf . .  At  the  same  time,  for  the  critical  states  the  corre¬ 
sponding  deflection  parameters  and  are  very  close 

or  coincide. 


For  a  twice  nonlinear  case  ( /#/,> ,  Fig.  1  shows  relationships 
Pi  (ut)  for  a  number  of  values  K 2  and  h/6  =0,005 3  and  Fig.  2  shows 
relationships  fn(u'6)  for  a  square  panel  with  =  40  and  k/S  =  o,002 
with  loading  areas  which  differ  with  respect  to  size. 

Here  r  -  parameter  of  the  total  load  per  area. 

/.  n 


obvious  algebraic  transformations  we  get: 


I  <’.  *  ’♦>  -  I  /,  (")  I  I  (*,  •?♦)  -  <7,  *,)  }  ?*  (0)  - 

/»Tn. /*./.) » » />;  <t,  •/*•?»>  »- 

»  /»' <3,  H>  ?*>  -  /»}  (»•  <-jA ).  ( 26 ) 

If  v/e  introduce  still  more  designations  to  certain  equations 

In  (*i  —  Ui '■**). 

In  <»,♦'*♦)  — 

(•»./*)  (i,  /I,  A)  1  it-Af/J,  ( 27 ) 

(jfflj’T*)  I  (?»'/VTi). 

1<VA)  -  Ivc,-M 

then  equation  (26)  is  rewritten  in  the  form  of 

/,••/* (0)  l-  •  ?*(0)  /': («,•/*)  I  f  A).  ( 28  ) 

By  performing  algebraic  calculations  for*  equation  (232)  an<3  the 
second  trio  of  vectors  $2,  72,  a2,  we  get  equation 

+  /.•:•?, (0)  («,•/,)  |. /,’(*,  + pi(Cyf+).  (29) 

Here  the  numbers  and  vectors  aOJ  SOJ  are  determined 

21 5  ^2’  2 ’  2  0 

*f  ^ 

by  the  same  formulas  as  numbers  and  vectors  a~,  b^,  C^, 

by,  of  course  substituting  index  1  in  these  formulas  for  the 
indicated  quantities  and  index  2  for  vectors  a,  I,  y. 
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From  Fig.  2  it  is  evident  that  the  localization  of  the 
load  on  a  smaller  area  leads  to  a  considerable  increase  in 
maximum  deflection  and  to  a  decrease  in  the  value  of  the  upper 
critical  load.  Thus,  with  a  load  distributed  along  the  central 
area  equalling  one  fourth  of  the  total  panel  area  (curve  pr  ) 
deflections  in  the  precritical  state  having  increased  almost 
three  times,  while  cracking  load  p*  has  decreased  by  approxi¬ 
mately  two  times,  as  compared  to  the  corresponding  values  for 
a  panel  uniformly  loaded  along  the  entire  surface  (curve  p,  ). 
With  a  load  of  pi  the  crack  virtually  disappears. 

Figure  3  shows  relationships  p,(:sc)  with  various  ratios 
h/&  for  a  square  panel  (  A  =i,  k-- =  Vo)  uniformly  loaded  along  the 

entire  surface.  It  is  evident  that  with  an  increase  of  this 
ratio  the  effect  of  physical  nonlinearity  on  deflections  and 
critical  loads  increases  significantly. 


Fig.  3. 
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